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全 摘要 本 文 利用 中 国 科 学 院 空间 科学 与 应 用 研究 中 心 的 瑞 利 激光 雷达 首次 观测 到 了 平流 层 地 形 重力 波 活 动 的 现 
AN， 结合 美国 国家 环境 预报 中 心 (NCEP) 的 全 球 预 报 系统 (GFS) 的 风 场 数据 分 析 了 该 地 形 重力 波 的 基本 参数 . 与 
惯性 重力 波 相 比较 ,地 形 重力 波 的 密度 扰动 没有 下 传 的 相位 ,在 同一 高 度 上 ,其 扰动 相位 保持 不 变 . 北京 空间 科学 
与 应 用 研究 中 心 瑞 利 激光 雷达 自 2012 年 开始 观测 实验 以 来 ,已 经 观测 到 多 起 地 形 重力 波 活动 事件 . 本 文 以 2013 
1 的 观测 数据 为 例 , 人 研究 北京 上 空 的 地 形 重力 波 活 动 , 并 结合 GFS 风 场 数据 分 析 了 北京 上 平流 层 地 形 
重力 波 的 波长 、 es 传播 速度 等 参量 . 通过 分 析 得 到 在 2013 年 11 月 11 日 北京 上 空 存在 一 列传 播 方向 为 北 偏 
西 52. 4 ,水平 波长 为 5.5 km, 平 均 垂 直 波 长 约 为 6.0 km 的 地 形 重力 波 . 

关键 词 ” 瑞 利 激光 和 雷达; 地 形 重力 波 ; 山地 波 
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Abstract Gravity wave activity in the upper stratosphere is investigated using density data 
retrieved from the Rayleigh lidar of National Space Science Center, Chinese Academy of Sciences. 
Combining the Rayleigh lidar data with the wind data of National Centers for Environmental 
Prediction (NCEP) Global Forecast System (GFS), we study a mountain wave observed on 
November 11, 2013. The parameters of this mountain wave, such as propagation direction and 
propagation speed, have been calculated. 

Gravity wave perturbations are extracted from 0.5 hX1 km density profiles. The relative 
density perturbations are expressed by oO (z) = (0(z)—p00 (2)) /pz), where plz) is the measured 
atmosphere density, and po (z) is the background density which is calculated by fitting the 
logarithmic form of whole night mean density with 4 order polynomial. The background wind 


data are achieved by applying a linear polynomial fitting to the NCEP-GFS wind data between 20 
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障碍 时 , 受 地 形 强迫 发 生 垂直 位 移 而 产生 的 波动 . 地 


地 球 物 理学 报 (Chinese J]. Geophys. ) 58 


to 48 km altitude. Using the data extracted from the complete density perturbations structure and 
the background wind data, we calculate the parameters of gravity waves observed on November 
11, 2013 by the gravity wave dispersion equation. 

The complete density perturbation structure shows an obvious phenomenon of mountain 
gravity wave activity. The wave phases at same altitude remain unchanged in the whole night. 
The perturbation structure shows that vertical wavelength is about 5.5 km but changes with 
altitude. A group of over-determined equations can be established by substituting the data 
extracted from the complete density perturbations structure and the background wind data into 
the gravity wave dispersion equation. And two groups of solutions are obtained by using the least 
squares method to solve these over-determined equations. The wind profiles in the direction of 
two sets of solutions have been analyzed. A critical layer (zero wind layer) which will prevent the 
upward propagation of mounting waves is found in the wind profile in the direction of 37. 9° (or 
217.9°). Finally, the gravity waves observed on November 11, 2013 propagate in the direction of 
52.4° from the north to the west, with a horizontal wavelength of 5.5 km. 

Compared with inertia waves, there is no downward-propagating or upward-propagating 
phase in the density perturbation structure. At the same altitude, the phase remains unchanged in 
the whole night. Such kind of gravity wave perturbation structures have been often observed in 
winter. Density data obtained by Rayleigh lidar and NCEP-GFS wind data at Beijing are used to 
analyze a mountain wave parameters observed on November 11, 2013. By analysis, we obtain 
terrain-generated gravity waves propagating in the direction of 52. 4° from the north to the west, 
with a horizontal wavelength of 5.5 km and average vertical wavelength of 6.0 km. 


Keywords Rayleigh lidar; Terrain-generated gravity waves; Mountain wave 


改进 算法 提高 信 噪 比 得 到 了 澳大利亚 和 非洲 的 好 
地 区 平流 层 地 形 重力 波 . 目前 ,国内 外 尚 没有 利用 
激光 雷达 对 地 形 重力 波 个 例 的 研究 . 本 文 从 地 形 
波 色散 关系 出 发 ,利用 中 国 


A 


引 二 


地 形 重 力 波 是 由 于 空气 流 经 山脉 .丘陵 等 地 形 
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瑞 利 
重力 


科学 院 空间 科学 与 应 用 


研究 中 心 的 瑞 利 激光 雷达 观测 的 数据 以 及 美国 


国家 


形 重力 波 对 平流 层 大 气 环流 、 热 力学 结构 以 及 化 学 
组 分 有 重要 的 影响 (Alexander et al. ，2011a). 目前 
国内 外 已 经 有 很 多 利用 卫星 探测 数据 对 平流 层 地 形 
重力 波 的 研究 工作 . Alexander 和 Teitelbaum (2007) 
通过 分 析 卫 星 探测 数据 ,发 现 2003 年 9 月 10 日 在 
南极 半岛 上 空 有 一 列 水 平 波长 为 300 km, 垂 直 波 长 
为 20 km 的 西北 向 传播 的 地 形 重力 波 . Jiang 和 Wu 
(2001) 利 用 UARS(Upper Atmosphere Research 
Satellite) 的 MLS(Microwave Limb Sounder) 数 据 
研究 了 安第斯 山脉 的 地 形 重力 波 , 曾 述 了 地 形 重力 
波 与 背景 风 场 的 相互 作用 ,指出 在 由 西风 主导 的 背 
景 风 场 下 ,地 形 重力 波 可 以 传播 到 50 km 的 高 度 . 
他 们 发 现 安第斯 山脉 处 的 地 形 重力 波 水 平 波 长 大 多 
为 110 km 和 400 km 左右 . 气流 流 经 澳大利亚 和 非 
洲 副 热带 地 区 的 山脉 形成 的 地 形 重力 波 波长 和 振幅 较 
小 ,不 易 被 卫星 探测 到 . Eckermann 和 Wu (2012) 通 过 


环境 预报 中 心 (National Centers For Environmental 
Prediction, NCEP) 的 全 球 预 报 系 统 (Global Forecast 
System,GFS) 的 风 场 数据 ,针对 北京 上 空 的 地 形 重 
力 波 个 例 进行 了 研究 ,并 分 析 计 算 了 该 重力 波 参数 . 


2 基本 理论 


地 形 重力 波 由 气流 流 经 地 形 障碍 (如 丘陵 ,山脉 
等 ) 而 形成 . 地 形 重 力 波 的 振幅 与 导致 其 产生 的 地 形 
障碍 成 比例 . 地 形 重力 波 相 对 于 地 面 是 静止 的 . 地 形 
重力 波 的 水 平 固 有 相 速 度 为 6 二 一 w，, 其 中 指 的 
是 与 地 形 相 垂直 的 背景 风 的 分 量 ; 波 的 水 平 传播 方 
向 与 水 平 风 速 方向 相反 ,如 果 风 速 随 高 度 变化 , 那 波 
速 也 会 随 高 度 变化 . 

地 形 波 的 固有 频率 为 : 


WwW = wo— kuo =— kuo. 
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根据 Fritts 和 Alexander(2003) 的 理论 推导 , 重 
力 波 的 色散 关系 为 : 
2 CN 1 
Wp A 
其 中 ,& 和 和 2 分 别 为 水 平 波长 和 垂直 波长 ，N 为 浮 
力 频 率 , f 为 科 里 奥 利 参 数 , 五 为 大 气 标高 . 


(2) 


在 mm 污 i ,太太 的 情况 下 ,色散 关系 可 
以 简化 为 : 
2 
m’ 一 a Ca 
Uo 


当 uo 二 0 时 ， 770 一 > SO， 此 时 ,垂直 波 速 为 零 . 我 
们 将 背景 风速 为 零 的 区 域 称 作 临界 层 ,在 临界 层 , 地 
形 重力 波 将 被 吸收 或 者 反射 

| | 二 || 时 , 加 为 虚数 ,重力 波 将 发 生 耗 散 


3 瑞 利 激光 雷达 实验 


中 国 科 学 院 空间 科学 与 应 用 研究 中 心 ( 简 称 空 
间 中 心 ) 于 2011 年 建立 瑞 利 激光 雷达 ,以 探测 30 一 
70 km 高 度 范围 内 大 气 的 温度 和 密度 . 该 瑞 利 激光 
雷达 的 参数 如 表 1. 

激光 光束 垂直 进入 大 气 , 在 30 km 以 上 的 高 度 
(一 般 认为 不 存在 气 溶胶 粒子 ) ,后 向 散射 光 主要 来 
自 于 大 气 分 子 的 瑞 利 散射 ,因此 后 向 散射 光子 数 正 
比 于 大 气 的 密度 . 大 气 的 相对 密度 值 可 以 由 光子 信 
号 直接 得 到 . 以 30 km 高 度 上 SABER 了 卫星 数据 的 
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表 1 空间 中 心 瑞 利 激光 雷达 参数 
Table 1 Parameters of Rayleigh lidar 


in National Space Science Center, CAS 


激光 波长 532 nm 
激光 功率 15 W 
激光 频率 30 Hz 
脉冲 宽度 7 ns 
脉冲 能 量 500 mJ 
望远镜 开口 面积 0. 4 m 


大 气 密度 为 参考 值 ,可 得 到 大 气 的 绝对 密度 (Wilson 
et al. ,1991; Marsh et al. ,1991). 大 气温 度 的 反 演 
在 理想 气体 以 及 静 力 平衡 的 假设 下 ,以 最 大 观测 高 
度 上 的 SABER 卫星 温度 数据 为 参考 值 , 利 用 大 气 
绝对 密度 得 到 大 气 的 温度 廓 线 . 温度 表达 式 为 
(Whiteway and Carswell ,1995 ) : 


ed | eadr 
plz) pp(z)Ra 
其 中 Rs 为 理想 气体 常数 . 
1 为 2013 年 11 月 11 日 观测 的 大 气温 度 廓 
线 , 累 积 观测 时 长 324 min, 距 离 分 辩 率 为 1 km. 其 
中 , 细 实 线 为 2013 年 11 月 12 日 北京 地 区 (40. 98"N， 
114. 10E) 的 卫星 温度 数据 . 


T(z) (4) 


4 重力 波 分 析 方 法 和 结 


4.1 利用 激光 雷达 数据 提取 重力 波 
空间 中 心 的 瑞 利 激光 雷达 自 建 立 开 始 , 已 有 49 


mm Measured temperature 
一 SABER temperature 


240 250 260 270 


Temperature/K 
图 1 2013 年 11 月 11 日 瑞 利 激光 雷达 观测 的 大 气温 度 以 及 SABER 温度 数据 
Hg.1 ， Atmospheric temperature measured by Rayleigh lidar and SABER temperature on November 11, 2013 
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个 夜晚 的 观测 数据 ,其 中 ,最 短 观 测 时 长 为 2 h, 最 长 
连续 观测 时 长 为 11 h, 累 积 观测 时 长 300 多 个 小 时 . 


重力 波 波 动 通过 瑞 利 雷 达 的 密度 数据 来 计算 . 
密度 扰动 的 表达 式 为 : p (zx) = 2 人 中 ,其 中 


plz) 为 观测 到 的 大 气 密度 , po (zx) 为 背景 密度 
(Hertzog et al. , 2001; Li et al. , 2010; Sivakumar 
et al. ,2006). 在 背景 密度 的 计算 中 ,考虑 到 密度 随 
高 度 指数 递减 ,因而 首先 对 整 晚 的 平均 密度 廊 线 取 
自然 对 数 再 进行 4 次 多 项 式 拟 合 ,进而 得 到 无 扰动 
的 背景 密度 曲线 (Gao et al. ,1998). 
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图 2 2013 年 11 月 11 日 9 
Fig.2 Half-hour relative density perturbation on November 11, 2013 


从 图 2 中 可 以 看 到 明显 的 重力 波 波动 ,并 且 在 
50 km 以 下 相同 高 度 上 的 相位 基本 保持 不 变 , 这 符 


合 地 形 重力 波 的 特征 . 该 重力 波 在 不 同 的 高 度 上 的 
垂直 波长 不 同 , 这 是 由 于 风速 随 高 度 变化 引起 的 . 通 
过 测量 波峰 与 波峰 之 间 以 及 波 谷 与 波 谷 之 间 的 距离 
可 以 得 到 不 同 高 度 上 的 垂直 波长 . 表 3 给 出 了 测量 
得 到 的 几 个 高 度 上 的 垂直 波长 . 

4.2 结合 GFS 风 场 数据 的 重力 波多 参量 分 析 方 法 


和 结果 
4.2.1 浮力 频率 
到 ve 2 lnb0 Es ao | 08 
频率 一 = ,其 
人 4 Bg 0 | eh 


c 一 VYR4sT 为 声速 ,7 为 定 压 热 容 与 定 容 热 容 之 比 . 
根据 浮力 频率 的 计算 公式 ,利用 当晚 的 观测 的 温度 
和 密度 数据 ,可 以 计算 得 到 浮力 频率 随 高 度 变化 廓 
线 , 如 图 3 所 示 . 
4.2.2 背景 风速 

背景 风速 会 影响 地 形 重力 波 的 垂直 波长 ,对 重 
力 波 的 传播 起 到 关键 作用 (Eckermann, 1995; 


以 2013 年 11 月 11 日 的 观测 数据 为 例 , 人 研究 北 
京 上 空 平流 层 的 重力 波 活 动 . 2013 年 11 月 11 日 的 
连续 观测 时 长 为 324 min, 用 于 计算 重力 波 扰动 的 
密度 数据 分 辩 率 为 0.5 hxl km. 表 2 给 出 了 四 个 
特征 高 度 上 半 小 时 积分 的 大 气 密 度 测 量 误差 . 
表 2 密度 测量 误差 


Table 2 Uncertainty of the density measurement 


高 度 (km) 30 40 50 60 
误差 (%) 0.0857 0.1842 0. 4235 0. 9722 


0.02 


0.00 


19:00 20:00 UT 


小 时 分 辩 率 的 相对 密度 扰动 


Altitude/km 


0.016 0.018 0.020 0.022 0.024 
Buoyancy frequency 


图 3 浮力 频率 


Fig. 3 Buoyancy frequency 


Collins et al. ，1996 ; 陈 操 ,2010). 在 分 析 计 算 地 形 
重力 波 参 数 时 ,需要 知道 当天 的 背景 风速 .我 们 利用 
NCEP GFS 分 析 资 料 给 出 的 北京 地 区 (40°N,116°E) 
2013 年 11 月 11 日 的 风速 来 计算 背景 风速 .GFS 分 
析 资 料 在 2013 年 11 月 11 日 有 0 时 .6 时、12 时 、18 
时 (UT 时 )4 个 时 刻 的 风 场 数据 . 由 于 我 们 的 观测 时 
间 为 16 时 到 21 时 ,因而 采用 18 ne 
20 一 48 km 高 度 范围 内 做 线性 拟 合 ,作为 背景 风 
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图 4 NCEP-GFS 风 场 数据 与 拟 合 的 风 场 数据 : 纬 向 风 
Fig.4 NCEP-GFS wind data and the fitted result:zonal wind 
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图 5 NCEP-GFS 风 场 数据 与 拟 合 的 风 场 数据 : 径 向 风 
Fig.5 NCEP-GFS wind data and the fitted result: 


meridional wind 


4 和 图 5 分 别 为 当天 20 一 48 km 高 度 范围 内 的 风 
速 拟 合 得 到 的 纬 向 风 和 径 向 风 结 果 . 
4.2.3 参数 计算 

根据 密度 的 扰动 图 ,可 以 大 致 测 出 不 同 高 度 上 
的 垂直 波长 . 我 们 将 不 同 高 度 上 的 各 个 参数 汇总 见 
表 3 所 示 . 


表 3 不 同 高 度 上 的 测量 参数 
Table 3 Parameters at different heights 


,wm 垂直 浮力 垂直 纬 向 径 向 ” 大气 

让 m) 波长 ”频率 波 数 风速 风速 。” 标高 
A(km) NCHz) mlm!) Un sr Vm*s 1) H(km) 
33.35 6.823 0.02259 9.209X10-4 8.97 人 6.39 
35,12- &00 0.02218 LL O47 X10°? 8. 56 he: 6.61 
39.02 4.385 0.02151 1.433X10™3 8. 25 2. 90 6. 94 
41.71 5.385 0.02094 1.167X10™3 2 yf 7.16 
44.50 5.769 0.02025 1.089X10™3 7.60 4.60 7.39 
48. 18 7.692 0.01928 8.168X10 14 7 5 人 和 
将 了 一 思 一 上 一 一 上 代入 到 色散 关系 式 

并 考虑 到 m? 六 Tr ,可 得 : 


kmiwu 十 有 220 — N) = mf. (5) 
在 表 中 任 选 两 个 高 度 上 的 数据 ,并 假设 波 的 传 
播 方向 角 为 90( 以 正 北 的 风 癌 为 0 ), 即 可 得 到 方程 组 : 


kiCmiw thkiu— Ni)= mf, (6) 
ki Cmzus t+ Riu — Ni) = mz f’, CL) 
ui = zolSin0 十 woicosO， (8) 
us 一 Wo Sing0 十 woz cOSO. (9) 


解 上 述 方程 组 可 得 到 重力 波 的 传播 方向 9 和 水 
平 波 数 & 以 及 水 平 波长 * = 2x/k. 对 于 每 两 个 不 同 
高 度 上 的 数据 ,可 以 得 到 两 组 解 . 表 4 给 出 了 传播 方 
向 0 和 水 平 波 长 4 的 平均 值 . 标 准 差 以 及 方差 . 
表 4 解 的 平均 值 \ 标 准 差 和 方差 


Table 4 Mean value, standard deviation and 


variance of the solutions 


平均 值 标准 差 方差 
01 40.5° 7.4° S54 12 
0. 122. 8° 6.6 43.9 
A1 5.7 km 1.2 km 1. 4 km? 
A2 5.6 km 0.7 km 0.5 km? 


我 们 还 可 以 采用 最 小 二 乘法 来 求解 上 述 超 定 


Wg 


方程 


表 3 中 ,任意 高 度 上 一 组 数据 代入 (5) 式 ,并 结 
合 (8)、(9) 式 可 以 得 到 一 个 包含 两 个 未 知 数 (&,0) 的 
非 线 性 方程 . 求解 这 样 一 个 非 线 性 方程 组 难度 较 大 . 
为 了 方便 计算 ,将 任意 两 个 高 度 上 的 数据 代入 (6) 一 
(9) 式 ,经 化 简 后 可 得 到 一 个 只 包含 一 个 未 知 数 (0) 
的 方程 . 由 6 个 不 同 高 度 上 的 数据 可 以 得 到 一 个 由 
15 个 方程 组 成 的 方程 组 . 然后 再 根据 最 小 二 乘 原 
理 , 求 使 得 该 方程 组 平方 和 最 小 的 解 0. 进而 再 求 得 
另 一 个 未 知 数 (8) 的 解 . 最 终 得 到 了 两 组 极 小 值 解 : 
0 =37.9°,k=0.00116 mi 一 5.4 km; 
0;=127.6°,k;=0.00114 m !,As=5.5 km. 
与 表 4 的 结果 相 比 较 , 用 最 小 二 乘法 得 到 的 结 
果 在 表 4 给 出 结果 的 误差 范围 以 内 . 两 种 方法 得 到 
的 结果 一 致 性 较 好 . 我 们 以 最 小 二 乘法 得 到 的 结果 
为 准 , 对 两 组 解 进 行 分 析 . 
为 了 检验 这 两 组 解 的 正确 性 ,下 面 分 别 给 出 在 
0 、0 方向 上 风速 随 高 度 的 变化 曲线 . 
由 图 6 可 知 ,0 方向 上 的 风速 在 低空 存在 零 风 
层 (临界 层 ) ,这 将 阻碍 地 形 重 力 波 向 上 传播 ,因此 ， 
观测 到 的 地 形 重 力 波 的 传播 方向 不 可 能 是 91. 而 在 
b 方向 上 ,从 地 面 到 48 km 高 度 范围 内 的 风速 都 为 
正 , 这 使 地 形 波 的 向 上 传播 到 达 上 平流 层 成 为 可 能 ， 
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— Windin oO 
--- Windinb, 
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0 10 20 30 40 50 
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图 6 0 和 方向 上 的 风速 
Fig.6 The wind profiles in O and 0 


考虑 到 水 平 波 数 k 的 方向 和 水 平 风速 方向 相反 , 且 
久 方向 上 风速 方向 为 正 , 所 以 & 的 方向 为 0 的 反方 
向 ( 即 307.6 ). 因 此 ,观测 到 的 地 形 重 力 波 的 传播 方 
向 为 307. 6〈 即 北 偏 西 52. 4 ) ,水 平 波长 为 5. 5 km. 


本 文 分 析 研 究 了 一 列 由 中 国 科学 院 空间 科学 与 
应 用 研究 中 心 瑞 利 激光 雷达 观测 到 的 北京 上 空 的 地 
形 重 力 波 , 这 是 国内 外 首次 由 瑞 利 激光 雷达 观测 到 
地 形 重力 波 的 案例 . 在 整 晚 的 密度 扰动 结构 图 中 可 
以 看 到 明显 的 重力 波 波 动 . 该 波动 在 不 同 高 度 上 具 
有 不 同 的 垂直 波长 , 且 相 位 在 同一 高 度 上 保持 不 变 . 
结合 NCEP-GFS 风 场 数据 ,我们 分 析 计算 出 了 地 形 
重力 波 基本 参数 ,如 波长 .传播 方向 等 参数 . 证 明 
2013 年 11 月 11 日 在 北京 上 空 存在 一 列传 播 方 向 
为 北 偏 西 52. 4° ,水 平 波 长 为 5.5 km, 平 均 垂直 波长 
约 为 6.0 km 的 地 形 重力 波 . 
致谢 ”背景 风 场 数据 获取 自 美国 国家 环境 预报 中 心 
(NCEP) 的 全 球 预 报 系统 (GFS) 网 站 . 
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